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ABSTRACT 

The observational data from some black hole candidates suggest the importance of electro- 
magnetic fields in the vicinity of a black hole. Highly magnetized disk accretion may play an 
importance rule, and large scale magnetic field may be formed above the disk surface. Then, we 
expect that the nature of the black hole spacetime would be revelled by magnetic phenomena near 
the black hole. We will start to investigate the motion of a charged particle which depends on 
the initial parameter setting in the black hole dipole magnetic field. Specially, we study the spin 
effects of a rotating black hole on the motion of the charged particle trapped in magnetic field 
lines. We make detailed analysis for the particle's trajectories by using the Poincare map method, 
and show the chaotic properties that depend on the black hole spin. We find that the dragging 
effects of the spacetime by a rotating black hole weaken the chaotic properties and generate 
regular trajectories for some sets of initial parameters, while the chaotic properties dominate on 
the trajectories for slowly rotating black hole cases. The dragging effects can generate the fourth 
adiabatic invariant on the particle motion approximately. 

Subject headings: magnetic field — black hole physics — relativity, chaos 

1. Introduction 

In the central region of AGNs, compact X-ray sources and GRBs, many black hole candidates are 
reported. The observational data from their active regions would include the informations about the gravi- 
tational and electromagnetic fields. However, there are no direct observational evidences for the existence of 
a black hole in the center of the compact active regions. So, it is very important how to find the informations 
on the curved spacetime from observational data, and how to interpret their nature. Although the inner 
part of a gas disk around a supermassive black hole emits non-thermal radiation, this may be explained by 
synchrotron emission from ultra-relativistic electrons in strong magnetic field regions around a black hole. In 
the central black hole system, however, the emission process is very little understood, although the consensus 
is in favor of emission due to plasma process. The emission from such particles would give us the information 
about the black hole spacetime and/or the electromagnetic field. 
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In a realistic model of a black hole-disk system, we should consider interactions between a plasma 
particle and a magnetic field in the curved spacetime. The plausible electromagnetic field around a black 
hole-disk system is still unkn own, although some authors have discussed this problem in the frame work of 
magnetohydrodynamics (e.g., Blandford fc Znaiekl Il977t ICamenzindl Il987t iNitta. Takahashi. fc Tomimatsu 
' 200l[ ). In general, the task for solving the structure of magnetic fields in a 
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black hole system is very hard. Then, to make a comprehensive analysis, we will give a vacuum external 
magnetic field around a rotating black hole, and investigate in detail the motion of a charged particle 
confined in the magnetic field . The sta t ionar y and axisymmetric vacuum ele ctror nagne ti c fields aroun d 
a black hole are discussed by 



PettersonI (|l975l) : IChitre fc Vishveshwaral (|l975h : (|2000l) : Ichoshl |2000l ): 



Tomimatsu fc Takahashil (|2001r ). We apply the black hole dipole magnetic field in Kerr geometry, which 
is a solution of the vacuum Maxwell's equations. Although the dipole magnetic field could correspond to 
the intrinsic field of a compact object like a neutron star, in the case of bl ack holes it has to results from 
current rings exterior to the event horizon but they can be very close to it ( Prasanna 19781 ). As known in 
the Earth's magnetosphere, the dipole magnetic field can trap charged particles within its magnetic bottle 
(i.e., the Van Allen belts). The particle gyrat es around the magnetic field line, and drifts in the toroidal 
direction (see, e.g., ICoedbloed fc Poedtsl l |200J). Furthermore, in the poloidal plane, the particle oscillates 



along the magnetic field line. In the black hole case, we can also expect similar situations. 

In a black hole magnetosphere, by solving equations of motion of a charged particle moving in the 
electromagnetic field numerically (see |Prasannalll978l . we can see the gyration around the magnetic field 
line, the drift in the toroidal direction and the multi-periodic motion in the poloidal plane. Such a orbit 
is very complicated in both poloidal and toroidal directions. Although a stationary and axisymmetric 
magnetosphere is assumed, there are two conserved quantities; that is, the energy and angular momentum 
of the charged particle. In addition to these, the rest mass of the particle is the third invariable quantity. 
However, the fourth invariable quantity related to the azimuthal motion of the particle does not exist when 
the electromagnetic field is considered around a black hole. Thus, the motion of a test charge in a black hole 
magn etosphere is not an integrable system. Then, the chaotic motion is expected (jNakamura fc Ishizuka 
19931 ). To examine the trajectory, we plot the Poincare map in the two-dimensiona l r-p^ plane, which shows 
inter sections of the particle's trajectories with the surface of section in phase space (jLichtenberg fc Liberman 
19921 ). If the plotted points form a closed curve, the motion is regular (not chaotic). This is because a regular 
trajectory moves on a torus in the phase space and the curve is a cross section of the torus. On the other 
hand, if the plotted points are distributed randomly, the motion is irregular (chaotic). From the distribution 
of the points in the Poincare map, we can judge whether or not the motion is chaotic. Then, we know chaotic 
behavior on the trajectory of a charged particle in a black hole magnetosphere. Furthermore, we find the 
black hole's spin effects on that chaotic motion. 

This study will be related to the source of radiation near a black hole. We will start this work as a basic 
step before considering the plasma as a magnetized fiuid. Although we find chaotic and/or regular motion 
in a black hole magnetosphere, the problem is how to find the phenomena in the observed spectrum. At this 
stage, however, we can not show the relation with the observational data in this work. The main purpose 
of this study is to understand the black hole spin effects (the dragging effects by a rotating black hole) on a 
test charge, which will be considered as the source plasma of high energy radiation in future. 

In §[21 we review the basic equations for the motion of a charged particle in a black hole magnetosphere, 
using the Hamilton- Jacobi formalism. To understand the nature of orbits of a charged particle in the black 
hole magnetosphere, we also consider in §[3] the effective potential including the four- vector potential of the 
magnetic field. We present the motions of a test charge in the black hole dipole magnetic field for various 
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values of black hole's spin parameters. In § |4] we make detailed analysis for its trajectories by using the 
Poincare map method, and show the chaotic properties that depend significantly on the black hole spin 
obviously. We see simultaneous presence of regular trajectories and regions of stochastic trajectories on the 
Poincare map. Furthermore, we find that due to the dragging effects of a rapidly rotating black hole, the 
fourth constant of motion approximately turns out to exit. Then, the motion of a test charge in a dipole 
magnetic field around a rotating black hole can behave like a integrable system. 



2. Equations for a Charged particle around a Black Hole 

We study the motion of a relativistic charged particle with the rest mass m and charge q in Kerr 
geometry. The background metric is written by the Boyer-Lindquist coordinates (t, r, 6, (p) with c = G = 1, 
where the signature of the metric is —2. The non-zero components of the contravariant metric g^'' are given 

by 



where E = + cos^ 6*, A = — 2Mr + a'^, A = (r^ + a^)^ — Aa^ sin^ 9. The external electromagnetic 
field is assumed to be stationary and axisymmetric. We assume that the electromagnetic fields would not 
disturb the background geometry of the spacetime. In fact the electromagnetic fields seen in astrophysical 
black hole candidates would be quite small compared with the gravitational field associated with the black 
hole. 



The Hamiltonian of the system can be defined as fe.g.. iMisner. Thorne fc Wheeleiill973f l 



H = ^crilT,, - qA,,)i7T, ^ qA,) (2) 

where tt^ is the canonical momentum and A^ = A^^ (r, 6) is the four- vector potential of the electromagnetic 
field. The four-momentum of a charged particle is 

P^^ — ^9^-{n,-qA,), (3) 

where A r/m is an afhne parameter and t is the proper time. The Lorentz- force equation is equivalent to 
Hamilton's equations written in terms of a;'* and tt^: 

dx" dH 

(4) 



dn^ dH 

d\ ^ dxf^ 



(5) 



In a stationary and axisymmetric black hole magnetosphere, the magnetic field is specified by a scalar 
function 'J of position called the stream function; that is, the poloidal magnetic field lines are contours of 
constant (The function ^! is proportional to the toroidal component of the vector potential, Af/,.) In this 
paper we set up a vacuum magnetosphere, so that no total poloidal current exists, and then the toroidal 
component of the magnetic field is zero. 

From the stationary and axial symmetry of both electromagnetic field and spacetime geometry, E = 
TTj — + qAt and L = —n^ = —{p^ + qA^f,) are constants of motion corresponding to the integrable 
coordinates t and (j) in eq. The third constant of motion is the particle's rest mass ra = {—g^^PfiPvY^^ ■ 
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In general, four constants of motion are needed to determine uniquely the orbit of a particle through four- 
dimensional spacetime. However, a test charge motion in a magnetic field around a black hole po sses onl; 
three obvious constants. So w e can expect chaotic behavior for its motions (Lichtcnberg fc Liberman.199 



Karas fc Vokrouhlick'vlll992f ). because such a system is non-integrable. Note that, without the magnetic 



field, the fourth constant of motion exists. This constant is known as Carter's constant of the motion 



Q^pI+ cos^ 



(6) 



which arises as a separation-of- va riables constant in the Hamilton-Jacobi derivation of equation of motion (see 
Misner. Thorne fc Wheeleiill973l ). That is, the motion of a test particle in Kerr spacetime without magnetic 
field is an integrable system. 

The stationary and axisymmetric electromagnetic fields around a rotating black hole in source-free re- 
gions are derived bv iPettersonI (119751 ). Here, we consider the solution that denotes dipole magnetic fields 
at distant regions, and call this solution "black hole dipole electromagnetic field" . In this black hole mag- 
netosphere, we expect that the dipole magnetic field configuration can trap a charged particle within their 
magnetic bottle, even if it extends close to the black hole. The four-vector potential of the black hole dipole 
has only two nonzero components, At and A^\ the dipole magnetic field and an induced quadrupole electric 
field in Kerr geometry are given by 



At 



Aa, = 



r{r ~ M) + (a^ - Mr) cos^ 



272I] 

- (r - A/ cos^ 
— 3/i sin^ d 



27 



472 E 



(r - M)a'^ cos^ d r{r^ + Mr + 2a^) 



{r(r^ - 2Mc? + aV) + Aa^ cos^ d\ — In 



27 \T — r- 



(7) 



(8) 



where 7 = (M^ 
axis. In the a — 



- 0^)^/2^ r± = M ± 7 and the dipole moment /i is taken to be anti-parallel to the rotation 
M limit, we obtain 



At 

Aa, 



-M/i[r sin^ 9 - 2{r - M) cos^ 9] 

2(r2 + A/2 cos2 9){r-MY ' 
-^sin^ 9[-2r'^ + M{r - M){r + M cof? t 
2{r^+M^ cos2 e){r-M f 



(9) 
(10) 



In the limit of a ^ M, \At\ and \A^\ diverge as (r — M)^'^, while in the case of a < M, \At\ and \A^\ 
diverge logarithmically at the event horizon. Although the source of the magnetic field should be outside 
the event horizon, the current loops could locate just outside the even t horizon (arbitrarily close to it) 
as the source of such dipole magnetic fields as considered in this work (|Prasannalll978) . In spite of this 
singularity, we will use this configuration for the study of the chaotic motion of a test charge around a black 
hole. The appearance of the singularity does not imply that the dipole field is invalid. It is valid in the 
regions considered outside the event horizon. In a realistic situation for the black hole magnetosphere as 
an astrophysica l model, to avoid the singularity at t he event horizon, the infinite sum of multi-pole fields 
IS necessary iTomimatsu fc Takahashill200l[ ). The plausible magnetic field configuration around a 

black hole should be investigated as a future work. 
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Orbits of a charged particle in a magnetosphere 



Now, we will discuss charged particle motions off the equatorial plane of the black hole dipole field. The 
concrete expressions of equations Q and ^ including electromagnetic field terms are very complicated, and 
are not particularly informative. However, by considering an effective potential in the poloidal plane, we 
can get a general picture of the orbits, although we need to integrate the equations of motion numerically 
to obtain the practical orbits. 

Here, we consider off-equatorial motion for a charged particle in the dipole magnetic field. By using the 
condition g^'^p^Pi, = with the constants of motion E and L, the equation of motion in the poloidal plane 
can be obtained. Then, without the kinetic terms of the poloidal motion {p^ — P^ — 0)i can define the 
effective potential ( Esteban fc Medina 1990l ) as 



E, 



nun ^ q_^^ 



5' 



L 

m 



-A 



5' 



P 



L 

m 



-I 1/2 



9 



(11) 



where i?min is the allowed minimum energy for a particle at the injection point. Although the distribution of 
the effective potential depends on the value of L/m and q/m of the charged particle and the electromagnetic 
field configuration At and A^j,, to study the chaotic behavior we are interested in the case that a charged 
particle is trapped in a bound orbit. When the energy of a charged particle is not so large (Kff < E/m < 1), 
the particle is bound in the orbit in the region of the potential well. That is, the particle has turning points, 
which correspond to the inner and outer envelope of the Larmor motion in the poloidal plane. 

The detailed motion i n the equa t orial plane (p^ = 0) in the Kerr background with the dipole magnetic 
field has been discussed by Prasanna ( 1980[ l. We will also see the projection onto the equatorial plane of the 
off-equatorial particle's motions. From the (^-component of the equation of motion, we obtain 



qAt) ~ g^^L + qA^) 



(12) 



For a non-rotating black hole case, the first term is negligible. When the signatures of parameters L and 
qfi are same, the particle does not gyrate because of p'^ ^ 0. To observe the gyration motion of a charged 
particle, where p*^ = is achieved periodically, we will choose L < and Qd = {q/m)fi > in the following 
numerical calculations. To see the off-equatorial motion of a charged particle, the initial value of p'' or p^ 
should be specified. In this paper, we inject the particle into the magnetosphere with the injection angle V'inj, 
which is the angle from the origin of the coordinate axes in the poloidal plane and is relate to the ratio of the 
initial p"^ and p^-values, where we choose < 4'ini < O.Stt for the calculations. The label "ini" indicate the 
quantity at the injection point. The equations are solved by using the 4th-order Runge-Kutta-Gill method. 

Figure [1] shows the hole's spin dependence on the effective potential Vcs{r,6) for the poloidal motion of 
a charged particle in the black hole dipole magnetic field. We see the valley-like structure in the lower-level 
potential region (named "potential valley"), which is across the northern and southern hemispheres nearly 
along a dipole magnetic field line, wherein the particle with lower energy can be trapped inside it. We also 
see a divide of the potential valley just on the equatorial plane. Figure shows the r- and ^-dependences 
of the local minimum of the effective potential in the r-direction, which is obtained from dVce (r, 0)/dO = 
and corresponds to the bottom along the valley, for various values of the spin parameter a. In this plot, we 
see the local minimum in the middle-latitude region in the ^^-direction. 

In the non-rotating black hole case (see Fig. [T^), the value of the local minimmn of the effective po- 
tential decreases from the equator to the event horizon along almost the dipole magnetic field line in both 
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( r/M) siQ e {r/M) sin 9 ( r/M } sin 9 

Fig. 1. — Effective potential with the black; hole dipole magnetic field (thin gray curves), (a) a — 0.0, (b) 
a = 0.3M and (c) a = M, where L/m = —7.0M and Qd = TG.OAf^. In these plots, the maximum level of 
Vcs is l.f and the interval is 0.02. 




Fig. 2. — Potential valley for various spin values, (a) (Kff)min vs. r and (b) (T4ff)min vs. 9 are plotted, 
where L/m = — 7.0M and Qd = TO.OM^. The circles in (a) indicate 6 = 7r/2, and the vertical line labeled 
by Ho. 04 shows the event horizon for a — 0.04Af . 



- 7- 



hemispheres, and has the minimum value at the event horizon. Thus, the potential valley is opened toward 
the black hole, while the width of the valley becomes narrower toward the event horizon. So, although a 
particle can be trapped within the potential valley for some time, the particle will fall into the black hole 
sooner or later. Figure [3] is an example of the orbit for a non-rotating black hole, where the projection 
on the equatorial plane (left panel) and the poloidal motion (right panel) are shown. The particle gyrates 
around the magnetic field line, and drifts in the toroidal direction. Furthermore, the particle oscillates in 
the poloidal plane along the magnetic field lines. In this dipole magnetic field case, the particle entering the 
regions of higher magnetic field strengths reflects back into the regions of smaller magnetic field strength. 
This is the so-called "mirror effects". Although the dipole field traps a particle between two mirrors, the 
particle falls into the black hole in the end, for almost non-rotating black hole cases. In Figure [3] we show 
the case of tpi^j = — O.IBtt, but it is easy to fall into the black hole sooner for the larger value of lAinj- 

Next, for a rotating black hole case (see Fig. [T]d), the centrifugal barrier by the hole's spin effects arises 
and the potential valley is disconnected from the event horizon. This is because the centrifugal barrier due 
to the hole's spin becomes higher with increasing the black hole spin. Thus, the dragging effects of spacetime 
enhance the mirror effects. Then, the particle can be trapped within the potential valley without falling 
onto the black hole. Furthermore, we see the "double-well potential" along the valley. Basically, the particle 
oscillates between the northern and southern hemispheres, but sometimes a particle with lower energy may 
be trapped in one hemisphere. Such a double- well potential may be related to the chaotic motion as discussed 
later. For the extreme rotating case shown in Fig. [TJ:, the spin effects are more effective. We see an almost 
single-well potential along the valley. Then, we can expect that the particle would oscillate between the 
northern and southern hemispheres periodically. The trajectory for the particle trapped in the potential 
valley will be discussed in ij5] again. 

Figures HHH show the orbits of a charged particle around a rotating black hole. Charged particles are 
trapped within the potential valley (i.e., a magnetic bottle) without falling onto a black hole. So, we can carry 
out long-time calculations for chaos studies. The particle oscillates quasi-periodically between the northern 
and southern hemisphere, although it is a very complicated orbit. In Figures [H [5l(T0P) and[6l(B0TT0M), 
the particle's orbit seems to be random, while in Figure [SKBOTTOM) and[SI^TOP), the orbit looks regular; 
that is, we can see some kind of periodicity despite its long-time motion. For rapidly rotating cases, we can 
see the regular orbits often in comparison with the slowly rotating cases. In general, some orbits show a 
regular trajectory, while a chaotic trajectory still remains. Thus, we can find spin effects on the motion of 
the charged particle. In the following, we analyze the detailed properties of these motions to clarify their 
dependence on the spin parameter. Such motions may be chaotic. So, we analyze the properties of the 
motions using a method with which we analyze chaotic behavior. 



4. Regular and Chaotic orbits 

The particle motions are common in the sense that they are the combinations of three types of motions 
(gyration, bouncing, and drifting). Here we analyze the property of such complicated motions using the 
Poincare map, which shows intersections of a trajectory with the surface of section in phase space. The 
motion of a point in phase space could be followed over hundreds of thousands of oscillation periods. The 
Poincare map is an useful tool to classify visually whether the motions are regular (non chaotic) or irregular 
(chaotic). To make the Poincare map, we adopt the equatorial plane {0 = 7r/2) as a Poincare surface and 
plot the point (r, p^) when the particle crosses the Poincare section with {p^ > 0). We obtain that, for a 
slowly rotating black hole case, a trajectory plots a lot of points on the map randomly, while for a rapidly 
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(/-/M)cos^ (r/M) sine 

Fig. 3. — The motion of a charged particle in the equatorial plane (LEFT) and in the poloidal plane 
(RIGHT). The particle finally falls into the event horizon from the polar region of the black hole. The 
parameters of the motion are given as a = 0, E/m = 0.920, Qd = 70M^, tpini = — O.IGtt, Xini — 10. 5A/ and 
^ini ~ n/2. In the poloidal plane, the dipole magnetic field lines (thin gray curves) and the effective potential 
(thin curves) are also plotted. 

rotating case there are some regular trajectories shown by tori. The different tori indicate the different initial 
injection angle ipini on the Poincare map. The appearance of regular trajectories suggests that the dragging 
effects of the black hole generate a nearly integrable system in spite of the existence of a magnetic field in 
Kerr spacetime (the details will be discussed in § [5|) . 

Figures [T UIOI show several typical examples of the Poincare map for various hole's spin and particle's 
energies. We find the spin effects on the trajectories in the Poincare map. For slowly rotating black hole 
cases, the region of random trajectories fill a finite portion of the energy surface in phase space; see Fig. [T] 
The intersections of a single random trajectory with the surface of section fill a finite area. For a rotating 
black hole cases, both types of trajectories are mixed; see Figs. [8] and [9l That is, the distribution on the 
Poincare map makes wide rings by random trajectories, and closed curves by regular trajectories. These 
differences on their trajectories depend on the initial ejection angle V'ini- Specially, for rapidly rotating black 
hole cases, the ratio of the regular orbits increases rather than mildly rotating black hole case; see Fig. [9] 
The hole's spin effects weaken the chaotic motion in the black hole magnetosphere. In fact, for the maximally 
rotating black hole case, for a wide range of ^ini-values the regular orbits are observed, although random 
trajectories also appear for larger energy particle motions; see Fig. [TUl 

Next, let us see the energy dependences on the Poincare map. For a slowly rotating black hole case (see 
Fig. [7|), the random nature of intersections on the Poincare map is almost independent of particle's energy. 
The outer boundary of cross sections plotted by the trajectory of pg"' = only shrinks with decreasing 
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Fig. 4. — An Example of the random trajectory of a test charge. LEFT panel shows the projections on 
the equatorial plane, and RIGHT panel shows the projections on the poloidal plane. The parameters of the 
motion are set as a = 0.3Af, E/m = 0.900, Vini = -0.207r, L/m = -7.0M, Qd = TO.OM^, nni = 10.5M and 
9ini - Tr/2. 



the particle's energy; it depe nds on the effective potential well specified by the values of Qd, L and a 



( Karas fc Vokrouhlick'v 1993) ■ We also see a few void regions on the map in some cases. On the other 



hand, we find the energy dependence of the Poincarc map for mildly and rapidly rotating black hole cases. 
In Figures [5^ and[Ht (where a — 0.6Af), most trajectories are chaotic, but in Figures [BJa and[HJi the regular 
trajectories appear on the map (less chaotic). Figure [H and (TU] show the cases of rapidly rotating black hole 
case. The regular trajectory is obtained when the elevation of the particle ejected is small (somewhat parallel 
to the equator) or large (somewhat perpendicular to the equator), while the random trajectory is observed 
when the elevation has an intermediate value of them. For example, in Figure [9^, when the initial value 
of is in the range of A-B or C-D, the regular trajectories are obtained, while the range of B~C shows 
random trajectories. 

Although the charged particle is trapped in the effective potential well, when the energy has an almost 
minimum value at the injection point; that is, E/m ~ K;ff (nni, 0ini), the intersections of the whole trajectories 
are bounded by outer and inner closed curves, which make a narrow layer; see Figs. [91 and[9ji. The inner 
boundary of the cross section is plotted for the orbit of = 0. Note that, in the case of a > 0.6Af, some 
orbits with their almost minimum energy show the random trajectories on the Poincare map, but the region 
of the random trajectories is restricted within a narrow torus-like belt. The width of this belt becomes 
narrow when the particle's energy has its minimum energy that is specified as the value of the effective 
potential at the particle's injected point; for a slowly rotating black hole case, such a belt-like distribution 
does not appear on the map. In summary, we see the tendency that random trajectories are independent of 
the initial injection angles and their energy in a slowly rotating black hole spacetime. For a rapidly rotating 
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(r/M) COS (r/M)smQ 

Fig. 5. — Examples of the trajectories of a test charge for a = 0.6AI. LEFT column shows the projections 
on the equatorial plane, and RIGHT column shows the projections on the poloidal plane. The parameters 
of the motion are set as (a,b) E/m = 0.895, 'i/'ini = — O.IOtt, (c,d) E/m — 0.890, V'ini = — O.SOtt, with 
L/m = -7.0M, Qd = TO.GM^, n^i = 10.5M and 0i„i ^ tt/2. 
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Fig. 6. — Examples of the trajectories of a test charge for a = 0.9il/. LEFT column shows the projections 
on the equatorial plane, and RIGHT column shows the projections on the poloidal plane. The parameters 
of the motion are set as (a,b) E/m = 0.885, i/'ini — — O.IOtt, (c,d) E/m = 0.885, 'i/'ini = — 0.227r, with 
L/m = -7.0M, Qd = TO.OM^, = 10.5A/ and e^in; = tt/2. 
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Fig. 7. — The Poincare map of a charged particle orbiting in the dipole magnetic field around a black hole 
of a = 0.3M. The intersections are shown by small dots on the map. The parameters for the motion are (a) 
E/m = 0.905, (b) E/m = 0.900, (c) E/m = 0.895 and (d) E/m = 0.892 with L/m = -7.0M, Qd = TO.OM^, 
rini = 10. 5A/ and 6'ini = it/2. 



- 13 - 



0.3 
0.2 
0.1 
0.0 
-0.1 
-0.2 

-0.3 



(a) 



1.0 9.0 10.0 11.0 12.0 13.0 

r / M 



0.3 r 

0.2 

0.1 

0.0 
-0.1 ■ 
-0.2 

-0.3 



(b) 



1 i 



8.0 9.0 10.0 11.0 12.0 13.0 

r/M 





8.0 9.0 10.0 11.0 12.0 13.0 

r/M 




8.0 9.0 10.0 11.0 12.0 13.0 

r/M 



Fig. 8. — The Poincare map of a charged particle orbiting in the dipole magnetic field around a black hole 
of a = 0.6A/. The parameters for the motion are (a) E/m = 0.905, (b) E/m = 0.900, (c) E/m = 0.895, (d) 
E/m = 0.890, (e) E/m = 0.885 and (f) E/m = 0.882 with L/m = -7.0M, Qa = TO.OAf^, n^i = 10. 5M and 
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Fig. 9. — The Poincare map of a charged particle orbiting in the dipole magnetic field around a black hole 
of a = 0.9M. The parameters for the motion are (a) E/m = 0.900, (b) E/m = 0.895, (c) E/m = 0.890, (d) 
E/m = 0.885, (e) E/m = 0.875 and (f) E/m = 0.872 with L/m = -7.0M, Qd = 70.01/2, ^ 10.5M and 
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Fig. 10. — The Poincare map of a charged particle orbiting in the dipole magnetic field around a black hole 
of a M. The parameters for the motion are (a) E/m = 0.900, (b) E/m = 0.890, (c) E/m = 0.880 and (d) 
E/m = 0.870 with L/m = -7.0M, Qd = 70.0Af^ nni = 10. 5M and ftni = n/2. 
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black hole case, however, we also see regular trajectories, which depend on the initial injected angle and also 
their energy. 



5. Discussion 

In § [21 we have mentioned that in a stationary and axisymmetric black hole magnetosphere there are 
only three constants of motion, vTf — E, t:^ — ~L and m, and as a consequence of the non-existence of the 
forth integral of motion, irg, the chaos would appear in the system. In spite of being a non-integrable system, 
however, we have found the regular trajectories on the Poincare map, although we have also seen the chaotic 
trajectories (for the different injection angles). In this section, we discuss the reason. 

As shown in Figures [JHH a particle in the black hole dipole magnetic field is trapped when a black hole 
rotates. It bounces back and forth between the northern and southern mirrors (also called magnetic bottle). 
Since the regular or random trajectories depend discontinuously on a choice of initial angle, their presence 
does not imply the existence of a global invariant of the system. However, if regular trajectories exist, they 
should represent some kind of invariants of the motion. These trajectories are conditionally periodic with 
angle variables. With this motion an adiabatic invariant would be associated, and that is rela ted to the 
longitudinal motion. Then, we define the action integral as (see, e.g., Goedbloed fc PoedtsI 2004 ) 



Je = ]<f V-Pep'^ dX , (13) 

where the integral is taken over one cycle of the oscillation in time and £ is the length of one cycle of the path. 
In general, the value of Jg is not a constant for each cycle of the oscillation of a test-charge in the black hole 
dipole magnetic field around a central object, so that we will see random trajectories (chaos). In fact, for a 
slowly rotating black hole, we see that most orbits with the various injection angles are chaotic, and their 
Je-values are not constants; they vary irregularly. However, we can also observe the regular trajectories, 
when the action integral of Jg is nearly constant for a long-time orbital motion. For a rotating black hole 
case, we can find the regular trajectory where the value of Jg becomes approximately constant, while it 
oscillates slightly in a few cycles. Note that the "Carter's constant Q" fluctuates in time; that is, the value 
of Q defined by equation ^ is not a constant. In this case, the motion shows a regular trajectory on the 
Poincare map. Thus, we find that the motion of a test charge in the magnetosphere around a rotating black 
hole can be nearly integrable. 

The region of random trajectories fill a finite portion of the surface of a section in phase space. FigurefTTk 
shows that the typical intersections of a single random trajectory with the surface of a section fill a finite 
area, where the value of Jg varies irregularly; see Fig. [TTb . The time series of 7r/2 — also shows an irregular 
feature in Fig. [TTb . In Figure [TTH . we see an annular layer randomly filled with intersections of a single 
trajectory lying between two invariant curves. We also see periodic trajectories that are composed by small 
regular trajectories inside these randomly filled regions (see, e.g.. Figs. [8ji and[9jl). Such a structure may 
be related to resona nces of the orbits that play a cr ucial role in the appearance of random motions in near- 
integrable systems ( Lichtenberg fc Liberman 1992h . In this paper, however, we do not discuss about the 



details. In this annular layer case, we also see that the value of Jg varies irregularly; see Fig. [TTb . On the 
other hand, in Figure [TTk . we see a generic trajectory that covers the surface of the torus. The motion along 
the dipole magnetic field line is almost periodic in the 6'-direction; see Fig. [TTl . The intersections of the 
trajectory with the surface of a section at a value of V'inj lie on a closed invariant curve and densely cover the 
curve over long periods of time. The value of Jg is nearly constant as shown in Figure [TTh . Figure [T2l shows 



-17- 





50 100 150 200 

cycle 




-15 



1000 2000 3000 

t/M 



1 1 ! 1 


(d) 






if 




■ 1 ^ 









8 9 10 11 12 13 

r/M 




50 lOD 150 200 

cycle 




-15 



1000 2000 3000 

t/M 





50 100 150 200 

cycle 




1000 2000 
t/M 



3000 



Fig. 11. — Various types of the Poincare map (LEFT column). Many points on the Poincare maps (a, 
d) are plotted by stochasticity (chaos), while that on the Poincare map (g) shows regular trajectory. The 
action integral Jq (CENTER column) and the time series of 7r/2 — 6 (RIGHT column), where the first 
several cycles arc shown, arc also presented. Small circles in the time series correspond to the points plotted 
on the Poincare map. The values of Jg are almost constant for (h) regular trajectories. The parameters 
of these trajectories are set as (a, b, c) a = 0.3M, E/m = 0.900, Vini = — O.IOtt, (d, e, f) a = 0.6M, 
E/m = 0.890, Vini = -O.IOtt, (g, h, i) a = 0.6M, E/m = 0.890, Vini = -O.SOtt. The remaining parameters 
are L/m = -7.0M, Qd = TO.OM^, ri„i = 10.5M and 6'ini = 7r/2. 
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Fig. 12. — Examples of the Poincare map (LEFT column) for a = 0.9M case. The regular trajectories are 
shown. The action integral Jg (CENTER column) and the time series of 7r/2 — 9 (RIGHT column), where 
the first several cycles are shown, are also presented. The parameters of these trajectories are set as (a, b, 
c) V'ini = — O.IOtt and (d, e, f) ipi^i = — 0.227r. The other parameters are E/m — 0.885, L/m — — 7.0M, 
Qd = lO.QhP, Tini = 10.5M and 6lini = 7r/2. 



the other characteristic trajectories in a rotating black hole spacetime; see also Fig. [51 In Figure [T^ TOP) . 
the crescent regular trajectory is presented, where the value of Je is approximately constant. On the other 
hand, in Figure [TSFBOTTOM). many islands of a trajectory distribute on a curve. This shows a regular 
trajectory. However, the value of Jg is not a constant, and it oscillates periodically. When we consider the 
average of Jg during dozens of cycles of oscillations, the value is almost constant over long periods of cycles. 
It seems that this feature is enough to guarantee the generation of the regular trajectory. 

The ratio of the regular trajectories on the Poincare map increases with the value of the hole's spin. 
Specially, for the maximally rotating case of a = M , almost trajectories become regular; see Fig. [TUl It 
seems that the fourth constant of motion is generated for the system with the black hole dipole magnetic 
field in the a ^ M limit. One may expect that, by using the dipole magnetic fields ^ and PH)) . the 
equation of motion can be separable with respect to the r and 6 coordinates, and the "modified Carter 
constant" including the magnetic field terms may be redefined. However, without the separable treatment 
on the basic equation, we find the constancy of Jg in this system numerically. Exactly speaking, for the 
regular trajectories appeared in this system the value of Jg is not necessary to be a constant; in fact it 
oscillates periodically. However, for a regular trajectory, we can regard the value of Jg (or the average on a 
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few cycles) as a constant approximately, comparing with that of chaotic trajectories. Even in the a — > M 
limit, this constancy can be broken depending on the injection angle. In fact, for larger E /m, we can also 
see an annular layer of randomly filled by trajectories. 

Our result suggests that the saddle in the double-well potential plays an important role for chaotic 
motions of a charged particle. Similar behavior has been investigated in Hamiltonian dynamical systems 



([Reichl fc Zheneill984( ). In the limit of a — > M, however, we see that the feature of the double well potential 
vanishes. Then, the chaotic behavior weakens, where many regular trajectories are observed. The property 
of chaos and/or regular trajectories in the dipole magnetic field around a rapidly rotating black hole is 
related to the hole's spin dependence on the eff'ective potential, which is generated by the combination 
of the electromagnetic force, the gravitational force and the centrifugal force. Although in this paper we 
consider the black hole dipole magnetic field, we could understand the basic properties of a charged particle 
by considering the distribution of the effective potential when we treat a motion of a test charge in any 
arbitrary electromagnetic field around a black hole. It is a future work to confirm the above suggestion by 
more detailed analysis. 



6. Concluding Remarks 

In this paper, we have discussed the off-equatorial motion of a test charge in the black hole dipole 
magnetosphere. The charged particle gyrates around a magnetic field line, drifts in the toroidal direction 
and oscillates between the northern and southern hemisphere, and then the orbits becomes very complicated. 
So, the numerical study by the Poincare map is effective. Then, we find the spin dependence of the trajectory 
on the Poincare map. That is, for a slowly rotating black hole case, the particle's trajectory shows chaos. 
More interestingly, we also find that, for a rapidly rotating black hole case, the chaotic behavior of the 
trajectories weakens and the fourth invariance of the motion, which is an adiabatic invariant related to the 
longitudinal motion, can be approximately generated. 

The chaotic motion in a black hole magnetosphere may be related to the origin of cosmic rays. The 
magnetic mirror concepts play a prominent role in space and plasma astrophysics. Examples are the Van 
Allen belt in the magnetosphere of the Earth, where high energy particles are trapped in a magnetic bottle 
and some kind of acceleration mechanisms are expected. In the case of a black hole magnetosphere, we may 
expect the formation of the magnetic bottle near the black hole. When such a situation is realized, the 
particles trapped in the "black hole Van Allen belt" would emit high-energy radiation, or some particles may 
accelerate to ultra-relativistic velocity by some electromagnetic interactions (ex., the Fermi acceleration by 
the perturbed magnetic bottle). Furthermore, some kinds of waves (i.e., Alfven waves and/or fast waves) 
in the magnetosphere will give influence to the motion of a particle moving in a magnetized plasma. Our 
simple model with the black hole dipole magnetic field presented here could be a preliminary step toward 
the understanding the motion of charged particles around a black hole. 

In this paper, we do not consider the interactions between charged particles, and the emission process 
from the chaotic/regular motion of a charged particle. However, we can expect that the spectrum emitted 
from such charged particles in periodic motions in the inhomogeneous magnetic field carries the informations 
on the black hole spin and the strength and/or distribution of the electromagnetic field. These studies should 
be considered in future works. 
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